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MODELING OF SUBSTITUTIONAL SITE PREFERENCE IN ORDERED 
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Abstract 

We investigate the site substitution scheme of specific alloying elements in ordered com- 
pounds and the dependence of site occupancy on compound stoichiometry, alloy concentration. 
This basic knowledge, and the interactions with other alloying additions are necessary in order to 
predict and understand the effect of various alloying schemes on the physical properties of a mate- 
rial, its response to various temperature treatments, and the resulting mechanical properties. Many 
theoretical methods can provide useful but limited insight in this area, since most techniques suf- 
fer from constraints in the type of elements and the crystallographic structures that can be mod- 
eled. With this in mind, the Bozzolo-Ferrante-Smith (BFS) method for alloys was designed to 
overcome these limitations, with the intent of providing an useful tool for the theoretical predic- 
tion of fundamental properties and structure of complex systems. After a brief description of the 
BFS method, its use for the determination of site substitution schemes for individual as well as 
collective alloying additions to intermetallic systems is described, including results for the con- 
centration dependence of the lattice parameter. Focusing on B2 NiAl, FeAl and CoAl alloys, the 
energetics of Si, Ti, V, Cr, Fe, Co, Ni, Cu, Zr, Nb, Mo, Ru, Hf, Ta and W alloying additions are 
surveyed. The effect of single additions as well as the result of two simultaneous additions, dis- 
cussing the interaction between additions and their influence on site preference schemes is consid- 
ered. Finally, the BFS analysis is extended to ternary LI 2 (Heusler phase) alloys. A comparison 
between experimental and theoretical results for the limited number of cases for which experi- 
mental data is available is also included. 
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Introduction 


During the past few years, a convergence of the computational sciences, condensed matter 
physics, chemistry and materials science has occurred in such a way as to provide valuable tools 
and expertise that could become of great utility in the design of novel, technologically useful 
materials. In the past, the development of new structural alloys and even small improvements to 
current alloys have mostly been performed through extensive experimental trial and error, which 
is both costly and time consuming. This approach is slowly changing, due to the concurrent 
effect of powerful computational techniques for different length scales, which provides basic 
knowledge needed at each step of new material development. 

Of the several computer-based tools now available to the materials scientist, atomistic simu- 
lations based on sound theoretical methods have the potential to become useful at the early 
stages of materials design, due to their ability to provide fundamental information at the atomic 
level. Atomistic studies are broadly based on two types of approaches, first principles methods 
and semiempirical techniques. These techniques are useful in separate but complementary 
areas. First principles methods provide detailed electronic structure information, while semiem- 
pirical techniques are good for large scale simulations for the determination of bulk structural 
properties. Ideally, first-principles approaches are best suited for providing the most accurate 
and consistent framework for any study at the atomic level. However, the complexity of most 
real life problems convoluted with the overwhelming computational requirements for even sim- 
ple systems have, to date, prevented first principles methods from becoming common and eco- 
nomical predictive tools beyond their current use. In fact, very few calculations exist that go 
beyond elementary binary alloys, and in some of these cases there are restrictions in scope and 
accuracy. This is gradually changing, and it is expected that future growth in computer power 
will correspondingly accelerate progress in this area. On the other hand, semiempirical methods 
based on quantum theory have enjoyed a great deal of attention in the last decade. This is mostly 
due to their computational simplicity, resulting in substantial progress in areas previously inac- 
cessible to theoretical treatment, such as alloy structure analysis and design. However, most of 
the methods developed in the past have been severely restricted in the type and complexity of 
the systems amenable to such studies. 
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The purpose of these semiempirical methods is to provide an efficient but accurate way to 
compute the total energy of arbitrary atomic systems in terms of the geometrical configuration of 
the atoms, using approximations for the interactions between atoms. When applicable, the results 
contribute a great deal to the general understanding of the system under investigation. However, 
in most cases, the existing techniques are restricted to a limited number of systems or they require 
specific (and therefore non-transferable) parameterizations or potentials. Naturally, it would be 
desirable to overcome these limitations on the predictive power and applicability of these meth- 
ods to increase their overall usefulness. 

Recently, a new semiempirical method was developed with the goal of avoiding these limita- 
tions so that it could be adopted as a materials design tool. The Bozzolo-Ferrante-Smith (BFS) ( 1 ) 
method for alloys, based on quantum perturbation theory, is particularly designed to deal with 
complex systems and geometries. It has no constraints in its formulation that would limit the type 
or number of elements under consideration or on the number or type of resulting phases. The two 
distinguishing features of BFS are its novel way of modeling the alloy formation process - by 
breaking the process into two coupled transformations - and its extremely simple mathematical 
formulation, requiring the solution of a single transcendental equation for each transformation. 
These features of the BFS method are particularly useful when dealing with multicomponent sys- 
tems, as will be shown in the present work dealing with the role of alloying additions to high tem- 
perature ordered intermetallics. 

It is now well recognized that alloying additions to structural intermetallics are essential for 
the optimization of physical, chemical and mechanical properties. However, the role of these 
additions in controlling properties is poorly understood due to the lack of detailed microstructural 
information. As an example, in a NiAl+Ti system, Ti additions in the order of 2.5 to 3.0 at. % have 
been shown to result in a 200 to 5000 fold reduction in creep rate as compared to that of the binary 
NiAl (2). However, from a microscopic point of view, this behavior is not clearly understood, thus 
limiting our ability to modify and further improve the alloy (3). 

In this work we take advantage of the versatility of the BFS method to tackle the problem of 
site substitutions of alloying additions to ordered intermetallics. Concentrating on transition metal 
rich B2 intermetallic alloys (NiAl, FeAl, CoAl), we studied the change in physical properties due 
to solid solution alloying additions, determined the substitutional site preference schemes, as well 
as interactions and changes due to multiple additives. To demonstrate the versatility of the BFS 
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method for this type of applications, we also analyzed the substitutional behavior of alloying addi- 
tions in a selected group of Ni 2 (Al,X,Y) Heusler alloys. 

The BFS Method 

Since its inception a few years ago, the BFS method has been applied to a variety of problems 
(1), ranging from bulk properties of solid solution fee and bcc binary alloys to more specific appli- 
cations like the energetics of bimetallic tip-sample interactions in an atomic force microscope 
and the effect of stoichiometry on the defect structure of NiAl (4) and FeAl alloys (5). The BFS 
method has also been used to deal with alloy surfaces (1). These previous studies provide a foun- 
dation for the work presented in this paper. 

An interesting consequence of the simple mathematical formulation of the BFS method is that 
simple expressions can be derived for predicting the composition dependence of some bulk alloy 
properties based solely on pure component properties (the BF rule) (6). As it will be shown later 
in this work, these tools will become particularly useful when dealing with multiple additions to 
binary or higher order alloys (e.g. Ni 2 AlTi). 

In what follows, we provide a brief description of the operational equations of BFS in order to 
introduce some concepts used later on. The reader is encouraged to seek further details in previ- 
ous papers ( 1 ,4-6) where a detailed presentation of the foundation of the method, its basis in per- 
turbation theory and a discussion of the approximations made are clearly explained. 

The BFS method provides a simple algorithm for the calculation of the energy of formation of 
an arbitrary alloy AH (the difference between the energy of the alloy and that of its individual con- 
stituents). In BFS, the energy of formation is written as the superposition of individual contribu- 
tions of all the atoms in the alloy, 


AH = 5>W.) = E e i. (1) 

i i 

where e t is the energy of an atom i in a pure elemental crystal and e\ is the energy of the same 
atom in the alloy being studied. For each atom, we partition the contribution e, to the energy of 
formation AH into two parts: a strain energy and a chemical energy contribution. The first contri- 
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bution takes into account the atomic positions of the neighboring atoms to atom i, regardless of 
their chemical identity. For its calculation, we use the actual geometrical distribution of the atoms 
surrounding atom i, computed as if all of its neighbors were of the same species as atom i. In this 
sense, the BFS strain energy differs from the commonly defined strain energy in that the actual 
chemical environment is replaced by that of a single-element crystal. Its calculation is then 
straightforward and even amenable to first-principles techniques. In our work, we use Equivalent 
Crystal Theory (ECT) (7) for its computation, due to its proven ability to provide accurate and 
computationally economical answers to most general situations. 

The chemical environment of atom i is evaluated in the computation of the BFS chemical 
energy term, where the surrounding atoms are forced to occupy equilibrium lattice sites corre- 
sponding to the reference atom i. Building on the concepts of ECT, BFS implements a straightfor- 
ward approach for the calculation of the chemical energy, properly parameterizing the interaction 
between dissimilar atoms. 

Thus defined, the BFS strain and chemical energy contributions take into account different 
factors, i.e. geometry and composition, computing them as isolated effects. A coupling function g 
restores the relationship between the two terms. The coupling function is defined in such away as 
to properly consider the asymptotic energy behavior where chemical effects are negligible for 
large separations between dissimilar atoms. Consequently, the contribution e, to the energy of for- 
mation AH from an individual atom i is then given by 

e i = e 5 + g£ C . (2) 


In what follows, we provide a brief description of the essential steps in the calculation of e, , 
which should be complemented with additional details described at length in Ref. 1. The strain 
energy contribution is obtained by solving the ECT perturbation equation (7) 


NR p x e aR 1 + MR p 2 e 


( a + l) R > 


= 5 >? 


■ e 


-(a + S(r ))rj 


(3) 
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where N and M are the number of nearest- and next-nearest neighbors respectively, and where p, l, 
a and X are ECT parameters that describe element i (see Ref. 7 for definitions and details), r 
denotes the actual distance between the reference atom and each of its neighbors and S(r) 
describes a screening function (7). The ECT parameters, together with the LMTO results, 
uniquely describe the physical properties of the element: p is related to the principal quantum 
number n ( p = 2n-2), a parameterizes the electron density in the overlap region between two 
atoms, 1 is a screening parameter that takes into account atoms located at distances greater than 
nearest-neighbor distances, and 1 is a scaling parameter that ensures that the dependence of the 
binding energy per atom as a function of lattice parameter follows satisfies the universal binding 
energy relationship of Rose et al. (8). The sum runs over nearest and next nearest neighbors. Eq. 
3 is solved for a s , the lattice parameter of the strain equivalent crystal, where the reference atom i 
has the same energy as it has in the geometrical environment of the alloy. Rj and R 2 denote the 
nearest- and next-nearest neighbor distances in this equivalent crystal and are therefore related to 
the unknown a s . Once the lattice parameter of the strain equivalent crystal a s is determined, the 
BFS strain energy contribution is computed using the universal binding energy relation of Rose et 
al. (8), which contains all the relevant information concerning a single-component system: 

e s = £ c ('l-(I+4)«^) < 4 > 


where Eq is the cohesive energy of atom i and where the scaled lattice parameter a $ is given by 


where q is the ratio between the equilibrium Wigner-Seitz radius and the equilibrium lattice 
parameter a e . The values of the equilibrium lattice parameter a e the cohesive energy E c and the 
bulk modulus Bq for the bcc phase of all the elements used in this work are displayed in Table 1 . 
These parameters are obtained from first principles, all-electron, density functional calculations of 
the elemental constituents in the symmetry of the alloy (i.e., bcc-Ni, bcc-Al, etc.). The particular 
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LMTO results 


Atom Lattice Cohesive Bulk 
Parameter Energy Modulus 

(A) (eV) (GPa) 


ECT parameters 





Ni 

2.752 

5.869 

249.2 

6 

3.067 

0.763 

0.2716 

Al 

3.192 

3.942 

77.3 

m 

1.8756 

1.038 

0.3695 

Cr 

2.837 

4.981 

286.1 

6 

HUB&a 

0.2300 

Ti 

3.213 

6.27 

121.0 

6 

2.821 1 

1.0476 

0.3728 

Cu 

2.8225 

4.438 

184.55 

m 

3.0492 

0.7615 

0.2710 

Co 

2.7347 

6.406 

285.81 

m 

3.4779 

0.7469 

0.2658 

Zr 

3.5525 

7.353 

93.77 

8 

3.4641 

1.2257 

0.4362 

Hf 

3.4942 

7.411 

113.02 

10 

4.3132 

1.1302 

0.4022 

Ta 

3.3091 

9.310 ' 

205.06 

10 

. 4.4230 

0.9664 

0.3439 

Nb 

3.3165 

7.844 

179.41 

8 

3.5099 

0.9472 

0.3371 

Si 

3.0701 

4.654 

99.99 

IB 

1.9773 

1.0158 

0.3615 

Fe 

2.7423 

6.071 

311.53 

6 

3.0939 

0.6955 

0.2475 

Mo 

3.1781 

7.195 

273.54 

8 

3.5050 

0.7506 

0.2671 

V 

2.9770 

7.512 

199.01 

6 

3.0050 

0.9290 

0.3306 

W 

3.1932 

9.953 

315.98 

m 

4.4395 

0.8194 

0.2916 

Ru 

3.0700 

7.664 

338.32 

8 

3.6116 

0.7086 

0.2522 


Table 1. Linear-muffin tin orbital (LMTO) results for the lattice parameter, cohesive energy and 
bulk modulus for the bcc phases of all elements listed. The last four columns display the result- 
ing Equivalent Crystal Theory (ECT) (7) parameters determined from the LMTO results, p is 
related to the principal quantum number n for the atomic species considered (p=2n-2 ), a param- 
eterizes the electron density in the overlap region between two neighboring atoms, A, is a screen- 
ing factor for atoms at distances greater than nearest-neighbor distance and / is a scaling length 
needed to fit the lattice parameter dependence of the energy of formation with the universal 
binding energy relationship of Rose et al. (8). 
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implementation used is the Linear-Muffin-Tin Orbitals (LMTO) method (9) in the Atomic Sphere 
Approximation. The ensuing ECT parameters p, l, X, and a e are also included in Table 1 . 

The BFS chemical energy is obtained by a similar procedure. As opposed to the strain energy 
term, the surrounding atoms retain their chemical identity, but are forced to be in equilibrium lat- 
tice sites of an equilibrium (otherwise pure) crystal i . The BFS equation for the chemical energy 
is given by 


NR P x ‘e a,Rl + MR P ‘e 




I 

k 


V 


,, Pi 

N ik r \ t e 


+ M lk r 2 e 




(6) 


where N ik and M lk are the number of nearest- and next-nearest neighbors of species k of atom i. 
The chemical environment surrounding atom i is reflected in the parameters a,*, given by 

<*ik = a i + A ki (7) 



Ni 

.. 

Al 

c^ 

Ti 

Cu 

Co 

Zr 

Hf 

Ni 



-.03588 

-.09062 

.01914 

-.05384 

-.12210 

-.12681 

Al 

.09160 


-.01524 

-.08649 

.05438 

.17971 

-.06367 

.25061 

Cr 

.22482 

-.01696 


-.06107 

.02449 

.09571 

-.11502 

-.08881 

Ti 

.49580 

.23399 

.07343 


.20565 

.67092 

-.05062 

.36712 

Cu 

-.01708 

-.04993 

-.01180 

-.07356 


.13897 

-.06767 

.13561 

Co 

MM 

-.06792 

-.02624 

-.09593 

-.05734 


-.06399 

-.13074 

Zr 

1.00755 

-.04861 

.40891 

-.01583 

-.04778 

-.05068 


.14942 

Hf 

.45749 

-.07141 

.01656 



.54726 

-.07487 



Table 2 BFS interaction parameters a xy (first subindex indicates the row, the second indicates the 
column) and a kx (in A’ 1 )for X-Y alloys (X, Y = Ni, Al, Cr, Ti, Cu, Co, Zr, Hf), determined by fit- 
ting the lattice parameter and energy of formation of the corresponding B2 compounds via 
LMTO. 
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where the BFS parameters A*, (a perturbation on the single-element ECT parameter a ( - used in Eq. 
3) describe the changes of the wave function in the overlap region between atoms i and k. Once 
Eq. 6 is solved for the equivalent chemical lattice parameter a^, the BFS chemical energy is com- 
puted using 


e 


c 



1 — (1 + a c 



( 8 ) 


where y, = 1 if a c > 0 and y- = -1 if a* c <0. The scaled chemical lattice parameter is given by 


_* ( Qt 

— C[' z * ' ^ 

h 

Finally, as mentioned above, the BFS chemical and strain contributions are linked by a cou- 
pling function g which describes the influence of the geometrical distribution of the surrounding 
atoms on the relevance of the chemical effects, and is given by 


where the scaled lattice parameter a* s is defined in Eq. (5). 

In this work we use the BFS interaction parameters, A ki , determined following the procedure 
outlined in Ref. 4. In order to provide parameters to the BFS method, we need to calculate the 
equilibrium properties of the elemental solid for the same symmetry as the compound to be stud- 
ied, since BFS is referenced to the ground state properties of the system in that symmetry. Once 
these parameters are computed, they remain the same for any other calculation involving any of 
these elements as related to the given crystal symmetry, requiring no further adjustment or 
replacement. A partial list of the most relevant parameters is included in Table 2. 
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Results and Discussion 


a. Single alloying additions to binary B2 compounds 

While many aluminide-based intermetallic compounds have high melting temperatures, 
excellent oxidation resistance and low density, they are limited in actual structural applications 
due to the lack of sufficient room temperature ductility and elevated-temperature creep strength. 
Ternary and higher order elemental additions can alter these basic properties and, therefore, a gen- 
eral understanding of the substitutional site preference schemes could be of great help in the alloy 
design process. In this section, we perform BFS analyses of the site substitution patterns of ter- 
nary and higher order additions to NiAl, FeAl and CoAl. 

The results quoted in this paper are based on BFS calculations of the energetics of a 72-atom 
cell in the AB B2 structure. This number of atoms is sufficient to study site substitutions of rele- 
vance to this work. A and B represent the two simple cubic sublattices of the B2 compound and X 
represents a ternary alloying addition. From the basic 72-atom computational cell we build a cata- 
logue of several atomic configurations, each describing a specific site substitutional scheme. We 
introduce a convenient notation for such configurations: X(A) denotes an atom X substituting for 
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Figure 1: (a) An X atom occupying an A1 site, (b) a Ni site with the Ni atom occupying the avail- 
able A1 site, at nearest-neighbor distance from each other , and (c) same as (b), but with the substi- 
tutional and antistructure atoms separated by a distance greater than the equilibrium nearest- 
neighbor distance. Ni and A1 atoms are denoted by black and grey circles, respectively. The alloy- 
ing addition is denoted with a solid square. 
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an atom A on the A sublattice, and X(A)B denotes the same case but with the displaced A atom 
occupying a site in the B sublattice. In this second case, we distinguish between the pair of defects 
X-A being nearest neighbors of each other (labeled NN in the figures) and the pair being sepa- 
rated by distances greater than that (labeled ‘far’ in the figures). Fig. 1 shows the simplest case, 
when just one atom (representing the alloying addition) occupies either an A1 site (X(A1), Fig. l.a) 
or a Ni (or Fe, Co in FeAl or CoAI alloys) site (X(Ni)Al), with the Ni atom forming an antistruc- 
ture defect by occupying an A1 site in a N^oA^ggjXj 39 alloy. Figs, l.b and l.c distinguish 
between the different relative locations of the substitutional atom and the antistructure atom. For 
simplicity and also because most alloying additions are made at the expense of Al, we will hold 
the transition metal composition constant at 50 at. % in most of the calculations unless otherwise 
noted. We could just as easily do the same calculations for substoichiometric transition metal 
compositions by taking into account vacancies in the calculations but space and actual interest in 
transition metal poor compositions preclude a discussion of those alloys in the present paper. 


? -°- 3 5 
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c3 

> 

O 

c 

o 

* w* 

03 


NisoCALig siXi 39) 
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& -°- 551 
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<3 -0.65| 

G 

m 


Fe Ru Co Cu Si Cr 


Zr Hf It 
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Figure 2: Energy spectrum for alloying additions to Ni 5 o(Al,X) 5 Q. The thick solid line corre- 
sponds to X(A1) defects (see Fig. l.a), whereas thin solid lines and dashed lines correspond to 
X(Ni)Al defects as shown in Figs, l.b and l.c, respectively. 
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The site preference of the alloying addition is determined by computing the energy of forma- 
tion of each cell and comparing the values obtained. However, determining the configuration with 
the lowest defect energy might not be sufficient for understanding the precise behavior of the 
additive. While it is true that the configuration with the lowest energy is the most likely to be 
found, it is important to consider the difference in energy between configurations describing dif- 
ferent types of substitution patterns. If the difference between X(A1) and X(Ni)Al is significant, it 
is to be expected that atoms X will always occupy A1 sites. Experimental difficulties in determin- 
ing the site preference of certain elements could be related, among other factors, to the size of this 
energy gap. Moreover, we will see later that the interplay between several simultaneous additions 
can drastically change the behavior observed when each addition is considered individually, mak- 
ing the analysis of higher energy configurations even more important. 

NiAl: In order to highlight the energy differences between various atomic arrangements, we 
present the results in the form of an ‘energy spectrum’. Fig. 2 shows a spectrum of the energies of 
formation of the configurations described in Fig. 1 (for one atom of the alloying addition in the 
72-atom cell), for several different elemental additions to NiAl. A wide range of behaviors is 
immediately apparent. In order to categorize the difference, we grouped the various elements con- 
sidered by the magnitude of the energy gap, G (defined as the difference in energy between the 
lowest energy level and its immediately higher energy level, relative to the energy of the lowest- 
lying state). In doing so, we find four distinct groups: (a) Fe and Ru, with a gap of less than 1%, a 
negligible energy difference that is well within the fluctuations in the parameters used for this cal- 
culation, (b) Co and Cu, showing a small but distinguishable gap, (c) Si,Cr, Zr, Hf with a well 
defined gap and (d) Ti, V, Nb, Mo, Ta and W with a gap that is over 20% of the energy of forma- 
tion of the defect. These last two groups are also characterized by the fact that all the additives 
considered mostly occupy A1 sites. 

In order to consider the possibility of changes in site preference with changes in concentra- 
tion, we performed similar calculations, also in a 72-atom cell, for a number of configurations, 
corresponding to a Ni^A^ 22 ^ 2.38 &U°y- It should be noted that even if this concentration is 
greater than the solubility limit for X in NiAl alloys, it is feasible to compute the energetics of X 
in solid solution in NiAl, in spite of the fact that a second phase might form at that concentration. 
The calculation is based on a predefined atomic distribution, which may or may not reflect the 
correct ground state structure for that composition. This issue has been addressed in our previous 
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work on the effect of Ti and Cr additions to NiAl (10). In this work, we disregard the possibility of 
precipitate formation for the sake of investigating general trends for this set of candidate alloying 
additions. 

Fig. 3 displays some of the possible configurations included in the catalogue of possible 
structures composed of two solute atoms in a 72-atom cell, showing only those sections of the cell 
affected by the presence of the alloying additions and the different resulting defect structures. Fig. 
4 shows the energy spectrum corresponding to this set of configurations. Because of the larger 


Y 50 A1 47.22 X 2.78 

Y = Ni, Fe, Co X= Ti, Cr, Cu, Hf, Ta, Nb, Zr, Si, Mo, V, W, Ru, Ni, Fe, Co 



2X(A1) nn 2X(Al) f 2X(Y)A1 nn 



2 X(Y)A1 nn 2X(Y)Al f [X(Al)+X(Y)Al f ] NN 



[X(A1)+X(Y)A1 nn ] nn [X(A1) +X(Y)Al NN ] f [X( Al)+X(Y) Al f ] f 

Figure 3: Configurations for two atoms in a 72-atom cell. A section of the cell containing the 
defect is shown. Y atoms (Ni, Fe, Co in NiAl, FeAl and CoAl, respectively) are indicated with 
solid black disks, A1 atoms with grey disks and the additive X with a solid square. 
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number of possible defects due to their relative location in the lattice, the spectrum now includes 
‘bands’ of states corresponding to the three main types of defects: the two solute atoms occupying 
A1 sites (thick lines), the two solute atoms occupying Ni sites with the displaced Ni atoms in A1 
sites (dashed lines), and a combination of both type of defects (thin solid lines). This energy spec- 
trum can be seen as an extension of the one shown in Fig. 2. A comparison with the single atom 
additions (lower concentration) provides information on the changes in the energy gaps due to the 
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Figure 4: Energy spectrum of the configurations shown in Fig. 3. For each element, the different 
energy states are grouped depending on the type of defect: two X atoms in A1 sites (solid thick 
line), one X atom in an A1 site and one in a Ni site (thin solid line) and both X atoms in Ni sites 
(dashed lines). In those cases where the X atom occupies a Ni site, an antistructure defect is also 
created. The different levels for the same type of defect correspond to the differences due to the 
relative location of the defects with respect to each other. 
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Ni 

sublattice 


No strong 
preference 


Ni5o(Al+X) 5 o 


A1 

sublattice 


Ti,Cr,Fe,Cu 



Zr,Nb,Mo,Hf,Ta,W 


Ti, V, Cr, Zr, Nb, 
Hf, Ta, W 



0. 1 eV difference 
between Fe(Al), 
Fe(Ni)Al 



Cr, Fe, V. 

Co 

(for x A1 < 42.8 at %) 



Ti, V, Cr, Zr, Nb, 
Hf, Ta, W ( both 
sites when unfilled 
with constituents) 



Fe, Ru, Cu (0.01 
eV difference 
between X(A1)/ 
X(Ni)Al) 


Experiment/Theory 




Fe, Ru, Cu, Si, Cr, 
Zr, Hf, Ti, V, Nb, 
Mo, Ta, W 


Theoretical model (11) 
Allaverdoba et al. 


Constitution diagrams (12) 


First-principles calculations (13) 
Fu and Zou 


Theoretical Model (14) Kao et al. 


Pseudo ground state model (15) 
Hosoda et al. 


ALCHEMI (16) Anderson et al. , 
X-ray (EXAFS)(17)Chartier et al. 
Atom probe field ion microscopy 
(18) Duncan et al. 


X-ray diffraction analysis (11) 
Allaverdoba et al. 


ALCHEMI (19)Field et al. , 
ALCHEMI (20)Cotton et al. 


ALCHEMI (2) Kitabjian et al. 


ALCHEMI (21), Munrow et al. 


BFS method (this work) 


Table 3 Comparison between experimental (light shaded cells) and theoretical results for site pref- 
erence of alloying additions in Ni(Al+X) alloys. 
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higher concentration of the alloying additions, and its effect on the site preference scheme. For 
example, a larger number of Ti atoms leaves the energy of the Ti(Al) level almost unchanged, 
while the Ti(Ni)Al levels (dashed lines in Fig. 2, dashed/dotted lines in Fig. 4) are clearly pushed 
toward higher energies, increasing the preference of Ti for A1 sites. Fig. 4 also shows that, for the 
atomic species considered and the range of concentrations studied, there is no change in the basic 
site preference scheme with concentration. 

A comparison of the BFS results with experiment and other types of theoretical calculations is 
provided in Table 3. In general, the BFS results are in excellent agreement with experiment and 
are consistent with other theoretical calculations. Unfortunately, few experimental results exist 
and most theoretical calculations have been performed for a rather small number of elements, 
limiting the comparison to the few cases available. Still, the available data is sufficient to put the 
BFS technique into proper perspective. 

There is agreement between all the theoretical models and the experimental results with BFS 
predictions for the site preference of Co, with the exception of the pseudo ground state model of 
Hosoda et al. (15), where a critical composition (x A , = 42.8 at. % Al) is determined for a change 
in site preference behavior. Above this critical concentration, Co is expected to occupy Ni sites. 
The site preference of Fe in NiAl has been the subject of several detailed studies due to the diffi- 
culty in establishing its behavior. In our work, Fe is found to have a very small value of G, with a 
slight preference for Al sites at zero temperature. This is in agreement with the detailed first prin- 
ciples study of Fu and Zou (13), who found a small value (0.1 eV) for the energy difference 
between the two types of defect. Several experimental studies (16-18) confirm the preference of 
Fe for Al sites. The BFS results for all the other alloying additions (Si, Cr, Zr, Hf, Ti, V, Nb, Mo, 
Ta and W) are in agreement with all the available experimental studies : Cr (1 1,19-20), Ti (2, 11), 
V (11, 21) andZr, Nb, Mo, Hf, Ta, W (11). 

For completeness, we summarize the numerical results for all the Ni(Al,X) alloy configura- 
tions in Table 4. For each entry, we list the energy of formation (in eV/atom) and the equilibrium 
lattice parameter of the resulting alloy. This value of the equilibrium lattice parameter is obtained 
by minimizing the energy of formation of every configuration with respect to isotropic expansions 
or compressions of the periodic cell. The changes in lattice parameter with elemental additions to 
Ni(Al,X) are shown in Fig. 5. The results obtained can be used to parameterize linear expressions 
for the dependence of the lattice parameter with concentration a(x). In order to free these expres- 



Table 4 Energies of formation (in eV/atom) and lattice parameter (in A) for the defect configurations shown in Figs. 1 and 3 for NiAl+X alloys. The 
shaded boxes denote the lowest energy state. The elements are grouped according the the value of the energy gap G (see text). 
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Figure 5: Lattice parameter as a function of concentration of several alloying additions X ( X = 
Si, Ti, V, Cr, Fe, Co, Cu, Zr, Nb, Mo, Ru, Hf, Ta,W) to NiAl. The values of |i (Eq. 1 1) are given in 
parenthesis. 


sions from their dependence on the end values (i.e. equilibrium lattice parameter of binary NiAl, 
a NiAi determined by first principles methods (4)), we normalize a(x) by a NiAb obtaining the gen- 
eral expression 

= 1 + jiix . (11) 

a NiAl 

where x is the concentration of the allowing addition. The values of h are included in Fig. 5. 

As mentioned before, the BFS method also allows for the derivation of simple expressions 
that provide an approximate value of n based only on the single-element properties (lattice 
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Energy of Formation (eV/atom) 


parameter, bulk modulus and cohesive energy, as listed in Table 1). The BF rule (6) for the lattice 
parameter of the alloy as a function of concentration of each element is 


Y* X i B i a i 

a(x) = -i ( 12 ) 


where a, and B t are the equilibrium lattice parameter and the bulk modulus, respectively, of ele- 
ment i in the symmetry of the alloy. In the particular case studied in this work, XAl 5 Q_ x Y x , with 
X = Ni, Fe, Co, Eq. (12) becomes 


a{x) = 


1 2 „ 

2 a x B x + 


(H* 


2 Al B Al + xa 2 v B 


Y U Y 


( 13 ) 


\ a X B X + Q - X^ a Al B Al + xa Y B Y 
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Figure 6: Energy spectrum for Fe (Al+X) alloys with alloying additions of X. See Fig. 2 for 
details concerning the type of defect structures modeled. . 
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from which the approximate value of u can be derived.These expressions do not take into 
account specific ordered structures, therefore leaving out the strong bonding characteristics of 
these B2 compounds. However, in the case of several alloying additions, it is obviously a practi- 
cal way of determining an approximation for the behavior of the lattice parameter with composi- 
tion. 

FeAl: Similar calculations were carried out for Fe(Al+X) alloys, following the same proce- 
dure used for NiAl. Fig. 6 displays the energy spectrum for a low concentration of the alloying 
addition X. The numerical results for all possible configurations, both for low and high concentra- 
tion, are listed in Table 5 As for NiAl, the different alloying additions are ordered according to the 
relative magnitude of the energy gap between the two main site occupancy choices. Ni and Ru 
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Table 5 Energies of formation (in eV/atom) and lattice parameter (in A) for the defect configurations shown in Figs. 1 and 3 for NiAl+X alloys. The 
shaded cells represent the lowest energy configuration. The alloying additions are grouped according to the value of the energy gap G, listed below 
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exhibit a negligible gap G <3 % while Co and Cu show a clear distinction between the two possi- 
ble levels ( 3 < G < 10 %). All the other elements display a clear preference for A1 sites, with Si, 
Cr, Zr and Hf having a distinct gap ( 10 < G < 30 %) and the rest (Ti, V, Nb, Mo, Ta, W) showing 
preference for A1 sites with excessively large gaps (G > 70 %). 

Completing the presentation of results for FeAl, Fig. 7 shows the change in lattice parameter 
with concentration for all the alloying additions considered in this work, also showing the corre- 
sponding slopes n of the linear dependence of the lattice parameter on the ternary addition with 
concentration. Again, the slopes of the linear expressions for a(x) were normalized to the LMTO- 
determined value for the lattice parameter of FeAl, 

1+pjc (14) 

a FeAl 

where x is the concentration of the alloying addition. The values of (I are included in Fig. 7. 

Unfortunately, very few theoretical or experimental results exist for Fe(Al+X) alloys. The 
preference of Ni and Co for Fe sites agree with the theoretical results of Allaverdova et al. 
(ll).They also predict the same site preference for Cu and Cr, although they admit that their pre- 
diction for Cr is questionable. The small energy gap found for Cu in this work could be an expla- 
nation of the ambiguity found in their work. That is not the case for Cr, where BFS predicts a 
large gap favoring Al sites, as found experimentally by synchrotron X-ray diffraction experiments 
by Khosla (22), confirmed by the work of Munroe et al (21), who established the preference of Cr 
for Al sites by means of ALCHEMI and X-ray methods. First principles calculations using Local 
Density Functional theory by Fu and Zow (13) found Ti and Cr preference for Al sites and Ni 
preference with Fe sites, in complete agreement with our predictions. There is also excellent 
agreement with the quasichemical model introduced by Kao et al. (14), which predicts V prefer- 
ence for Al sites and Co and Ni preference for Fe sites. The authors point out that the predictions 
for other transition metal additions have large uncertainties and should not be used as basis of 
comparison. 

One of the most comprehensive studies of site distribution of transition metals is the recent 
work of Anderson (23), which includes an ALCHEMI analysis of Ti, V, Cr, Mn, Co, Ni and Cu in 
FeAl. The BFS results for Co, Ni and Cu exhibit a very small energy gap G. The slight preference 


111 T 


22 



of Co and Ni for Fe sites is in agreement with the ALCHEMI predictions, but Cu results are not. 
In these situations (i.e. small G), as noted before for Fe in NiAl, just the inclusion of entropic 
effects at finite temperature may or may not be sufficient to alter the preference scheme. 

The most marked preference for A1 sites is that of Ti, in agreement with our work, with 85 % 
of the alloying element occupying A1 sites. The residual Fe site occupancy is explained by Ander- 
son in terms of site equilibration kinetics, as diffusion of substitutional defects is required. Unfor- 
tunately, we have not included the possibility of vacancies in our calculation, as NiAl - the basis 
of our study- does not have structural vacancies for Ni-rich alloys. A more general catalogue, 
including vacancies, has been defined and tested in our previous work on binary FeAl (5). The 
defect structure (substitutions for Fe rich, triple defects for FeAl and Al-rich FeAl) was correctly 
determined, in agreement with experiment and other theoretical approaches. 
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Figure 8: Energy spectrum for Co(Al+X) alloys with alloying additions of X. See Fig. 2 for 
details on the type of defect structures modeled. 
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Table 6 Energies of formation (in eV/atom) and lattice parameter (in A) for the defect configurations shown in Figs. 1 and 3 for Co(Al+X) alloys. The 
shaded cells denote the lowest energy state. The elements are grouped according to the value of the energy gap G (see text). 
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Figure 9: Lattice parameter as a function of concentration of several alloying additions X to 
CoAl. The values of |i (Eq. 14) are given in parenthesis. 


CoAl: Having found agreement with available experimental and theoretical evidence for both 
NiAl and FeAl raises confidence in the BFS method and the parameterization used in this work. 
The parameters determined for the individual elements, as well as the BFS parameters a ab calcu- 
lated for all pairs A-B are unique and transferable, making it trivial to use them for the calculation 
of any other B2 alloy formed by the elements considered. It is therefore reasonable to extend the 
application of this methodology to CoAl, a system for which we could find no experimental data 
on the behavior of alloying additions in this compound in the current literature. Following the for- 
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mat of our previous examples, we display the BFS predictions for the site occupancy in Co(Al+X) 
alloys by showing the energy spectrum for single-atom substitutions in Fig. 8, and the complete 
set of results for the energy of formation and equilibrium lattice parameter as a function of com- 
position in Table 6. The elements are grouped by the value of their energy gap G: Fe, Ni (G < 2 %) 
; Cu, and Ru (2 < G <5 %); Si, Cr, Zr and Hf (5 <G<10 %) and Ti, V, Nb, Mo, Ta and W (G >20 
%) In addition, the change in lattice parameter is displayed in Fig. 9 (see Eq. 11). 

To conclude the analysis of ternary additions to XA1 (X= Ni, Fe, Co), we summarize the site 
preference scheme in Fig. 10, where we show, for each intermetallic compound, the site prefer- 
ence distribution. It should be noted that while helpful in obtaining a quick feel for the behavior of 
each element, there is no explicit consideration of the strong ordering tendencies that characterize 
these compounds. It is therefore to be expected that the BF rule predictions of the absolute values 
of the lattice parameter might not be as reliable as those obtained from full BFS calculations. As 
with the lattice parameter in Eq. 12, the bulk modulus, for arbitrary number or additions and con- 
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centrations, can be computed for any of the alloys studied or for those that can be constructed 
from the elements parameterized in this work, by using the following expression 



i 


where the single element values a, are listed in Table 1 . The values of 5, can be obtained from the 
definition of the scaling length /, 



\2%qB i a e 


(16) 


b. Double alloying additions to B2 compounds 

Beyond examining all possible B2 alloys that can be composed of the 16 elements considered 
in this work (i.e. CuAl, etc.), the BFS method allows for a straightforward generalization to higher 
order additions, without any new determination of parameters and without any additional com- 
plexity in the calculations or the computer resources needed. A review of the description of the 
BFS method should convince the reader that no particular provision should be taken when 
increasing the number of elements in the system. The only additional effort consists of generating 
a larger ‘catalogue’ of possible atomic configurations than that used for the simple ternary case. 
We therefore constructed a sufficiently large catalogue of atomic distributions contemplating all 
the possible atomic arrangements and the ensuing defects for two or more simultaneous alloying 
additions to binary B2 intermetallic compounds. One of the many goals of such effort is to inves- 
tigate the changes, if any, in site preference of individual atomic species in the presence of other 
alloying additions. 

As can be seen in Fig. 1, some types of atoms might be susceptible to changing their substitu- 
tional behavior in the presence of other additions due to the small energy gap between different 
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configurations. Ti, for example, has a definite preference for A1 sites in NiAl, a fact clearly 
proven experimentally and by any other theoretical calculation (2, 11, 12). The large energy gap 
between Ti(Al) and Ti(Ni)Al substitutions clearly support this fact. It is to be expected that the 
addition of a fourth element would be, in general, irrelevant in changing the behavior of Ti in 
NiAl. This can be seen in Fig. 11, which shows the energy spectrum of Ti in NiAl, and the 
changes it undergoes in the presence of a fourth element. A detailed analysis of this figure indi- 
cates that, in general, the Ti(Al) choice is not affected by any of the additions. The lower part of 
the spectrum is completely dominated by Ti(Al) defects (thick solid and dashed lines). Con- 
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Figure 11: Energy spectrum of two simultaneous alloying additions to Ni(Al+X+Y) alloys. The 
first column shows the ternary case, Ni(Al,Ti) The subsequent columns describe the energetics 
of quaternary alloys with the defect structure indicated in the inset. 
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versely, the presence of Ti does not seem to affect the individual preference of the other elements. 
Fe, for example, has a minuscule energy gap when acting alone, and it continues to do so in the 
presence of Ti: Ti(Al) +Fe(Al) and Ti(Al)+Fe(Ni)Al defects are almost equivalent in energy. 

The strong tendency of Ta, Mo, V and W for A1 sites is emphasized in the presence of Ti, 
making Ti(Ni)Al substitutions slightly more favorable than Ti(Al), although they are all very high 
in energy to have any significance. The numerical results are displayed in Table 7, which in addi- 
tion to the results plotted in Fig. 1 1 for Ni(Al+X) alloys, includes results for two-element addi- 
tions to stoichiometric NiAl. 

Co, on the other hand, benefits from the presence of Fe increasing its likelihood for choosing 
Ni sites with respect to its isolated behavior. The same effect is seen in Ru, but in favor of A1 sites. 

We could continue this analysis for all the possible combinations of alloying additions consid- 
ered in this work, as well as extending it to include 3, 4 or more simultaneous additions, a calcula- 
tion that would only entail the generation of the necessary catalogue of configurations. Such 
massive presentation of results is naturally beyond the limitations of this paper, and it will be pub- 
lished elesewhere. The main goal of this paper is to simply introduce the BFS method for this type 
of application with an appropriate illustration of typical results 

c. Heusler structures 

We conclude our analysis by considering single alloying additions to ternary Heusler phases. 
Once again, the complexity of the calculation depends on our ability to define the proper cata- 
logue of configurations which will adequately describe the system at hand. We continue to use the 
same set of parameters that were used for the B2 compounds. While numerous Heusler phases 
exist, there is very little known about these alloys in terms of their physical properties, except for 
the lattice parameter of some alloys. No information on site preference of alloying additions 
seems to be available. Fig. 12 shows the most important site substitution schemes consistent with 
the structure of the Heusler phase (A atoms in one cubic sublattice, and B and C atoms in the 
other sublattice, in alternating comers of the cube). Table 8 lists the corresponding energies of 
formation and equilibrium lattice parameter for each one of these atomic arrangements, from 
which the site preference scheme can be extracted for a base Heusler compound of N^AlTi. 
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Ti(Ni)Al+X(Ni) 

-.3971 

2.8849 

-.4103 

2.8830 

-.3242 

2.8941 

-.4462 

2.8686 

-.4577 

2.8657 

-.4380 

2.8687 

-.4181 

2.8712 

-.4359 

2.8724 

Ti(Ni)Al 

+X(Al)Ni 

-.3971 

2.8849 

-.3669 

2.8893 

-.3234 

2.8939 

-.4081 

2.8721 

-.4027 

2.8711 

-.4021 
2.871 1 

-.3886 

2.8745 

-.3870 

2.8773 


Table 7 Energies of formation and equilibrium lattice parameter for Ni 5 o(Al,Ti,X) 5 o and for Ni 5 o_ 
y/ 2 Al5o_ y /2(Ti,X)y alloys (where the Ni and A 1 ratio is kept constant at 1:1), for X = Hf, Zr, Ta, Cu, 
Co, Fe, Cr and Si. The shaded cells indicate the lowest energy state. The shaded cells indicate the 
lowest energy state for each element. In all cases it is assumed that the Ti and X atoms are located 
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(Ni,X) 2 AlTi Ni 2 (Al,X)Ti Ni 2 Al(Ti,X) 



Figure 12: Atomic distributions of X additions to a Ni 2 AlTi alloy. Ni-poor, Al-poor and Ti-poor 
alloys are indicated in each column. Ni, A1 and Ti atoms are indicated with solid black disks, 
solid grey disks and open circles, respectively. The additional element is denoted with the sym- 
bol X. 
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Ni 2 AlTi+X 

Hf 

Zr 

Ta 

Cu 

Co 

Fe 

Cr 


(Ni^) 2 AlTi 

X(Ni) 


- 0.3599 

2.9398 

- 0.2274 

2.9540 

g 1 ; y, 

|g?fpl§ 


IP? 

li:f ii : A . 

... . .j. 



X(AI)Ni 

- 0.3856 

2.9391 

•0.4102 

2.9373 

-0.3757 

2.9392 

-0.4314 

2.9226 

- 0.4130 

2.9230 

- 0.4144 

2.9227 

-0.4135 

2.9243 

X(Ti)Ni 

- 0.2782 

2.9488 

- 0.3043 

2.9467 

- 0.2845 

2.9471 

- 0.3238 

2.9322 

- 0.3058 

2.9324 


- 0.3016 

2.9345 

X(Ti,Al)Ni 

- 0.3810 

2.9396 

w* 

ipp- , 

-Ip- ; ! ||i = 

- 0.4268 

2.9231 

- 0.4074 

2.9235 

- 0.4089 

2.9232 

-0.4095 

2.9246 

X(Al,H)Ni 

- 0.2782 

2.9488 

- 0.3000 

2.9472 

- 0.2731 

2.9483 

m 

- 0.3067 

2.9324 

- 0.3084 

2.9320 

- 0.3013 

2.9345 


Ni 2 (Al,X)Ti 

X(A1) 

lljM’feri’S 

-0.4476 

2.9298 


ilSfl 

tMjifci • ii. : ! i : 

niiLL a- 

- 0.4366 

2.9164 

-0.4429 

2.9154 

| lj|T ]f ' V V i Ur f «Ibi 

i >: 

.. 

X(T1)A1 

-0.4466 

2.9293 


E. 

;; • 15*: :f 1 ^ ■ 

-0.4525 

2.9163 

- 0.4319 

2.9169 

-0.4384 ! 

2.9159 

-0.4348 

2.9178 

X(Ti,Ni)Al 

- 0.2842 

2.9426 

- 0.3104 

2.9405 

- 0.2897 
2.941 1 

- 0.3304 

2.9260 

- 0.3124 

2.9263 

- 0.3143 

2.9259 

- 0.3081 

2.9284 

X(Ni)Al 

- 0.4153 

2.9294 

- 0.3522 

2.9351 

- 0.2005 

2.9512 

m 



- 0.4101 

2.9194 

X(Ni,T1)Al 

- 0.4107 

2.9298 

- 0.3508 

2.9352 

- 0.2003 

2.9513 

- 0.4310 

2.9180 

-0.4596 

2.9138 

•0.4480 

2.9154 

- 0.4063 

2.9198 


Ni 2 Al(Ti,X) 


X(Ti) 

X(Al)Ti 



sipp 

•• JLwpsPp 

SifcS' 1 <•**. V’U"”; i i 

W&JXZviyrl 

- 0.4392 

2.9133 

- 0.4435 

2.9129 

-0.4457 

2.9123 


X(Al,Ni)Ti 

- 0.3779 

2.9317 

- 0.4040 

2.9297 

- 0.3720 

2.9314 

- 0.4245 

2.9152 

- 0.4058 

2.9156 

2.9153 

- 0.4073 

2.9167 

X(Ni)Ti 

- 0.4162 

2.9266 

- 0.3709 

2.9306 

- 0.2282 

2.9460 

- 0.4394 

2.9144 

-0.4658 

2.9104 

-0.4520 

2.9122 

- 0.4127 

2.9164 

X(Ni,Al)Ti 

- 0.421 1 
2.9261 

- 0.3745 

2.9303 

- 0.2281 

2.9460 

- 0.4435 

2.9140 


f !*1 1 1 .f 1 ' " 1 . ,J *• ^ 1 

- 0.4154 

2.9162 


Table 8 Hf, Zr, Ta, Cu, Co, Fe and Cr additions to Ni 2 AlTi alloys. Energies of formation (in eV/ 


atom) and equilibrium lattice parameter for the atomic configurations displayed in Fig. 12. The 
dark shaded cells indicate the lowest energy case. The light shaded cells indicate those states that 
are both close in energy and lattice parameter to the lowest energy state. 
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X(Ni) 


X(Al)Ni 


X(Ti)Ni 


X(Ti,Al)Ni 


X(AI,Ti)Ni 


X(A1) 


X(TI)A1 


X(Ti,Ni)Al 


X(Ni)Al 


X(Ni,Ti)Al 


X(Ti) 


X(Al)Ti 


X(Al,Ni)Ti 


X(Ni)Ti 


X(Ni,Al)Ti 


(Ni,X) 2 AlTi 


- 0.2340 

2.9524 


- 0.2580 

2.9482 



- 0.2891 

2.9457 


- 0.2924 

2.9449 


Ni 2 (Al,X)Ti 



- 0.3036 

2.9388 


- 0.2260 

2.9476 


- 0.2254 

2.9477 


- 0.4354 

2.9253 


- 0.2434 

2.9449 


- 0.2417 

2.9451 


- 0.3776 

2.9252 


Ni 2 Al(Ti,X) 


- 0.2227 

2.9481 



- 0.2511 

2.9426 


- 0.2509 

2.9426 


- 0.26000 

2.9407 


- 0.3851 

2.9221 


Table 8 (continued) Nb, Mo, V, W, Ru and Si additions to Ni 2 AlTi alloys. Energies of formation 
(in eV/atom) and equilibrium lattice parameter for the atomic configurations displayed in Fig. 12. 
The shaded cells indicate the lowest energy case. The dotted cells indicate those states that are 
both close in energy and lattice parameter to the ground state ones. 
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Energy of Formation (eV/atom) 


(Ni,X) 2 AlTi 

Ni 2 (Al,X)Ti 

Ni 2 Al(Ti,X) 


-0.40 
-0.41 — 
-0.42 — 
-0.43 — 
-0.44 — 
-0.45 — 
-0.46 I_ 
-0.47 — 


Hf(Ni) 


Hf(Ti)Al 

Hf(Al) 

Hf(Ti) 

Hf(Al)Ti 



Zr(Al) 

ZrfTOAl 

Zr(Al)Ti 

' ZrOi) ' 




CuCADNi 

Cu(Ni) 


Cu(Ti)Al 

Cu(Al) 

_ Cu(Ti)_ 
_Cu(A0Ti 


Fe(Ti)Al 


Fe(Ni)Al 

Fe(Ti) 

Fe7Ni.fT)A l 

_Fe£AlJTi_ 

je(Ni)TiI 

Fe(Ni)Al 


Co(Ni,Ti)A l FeTNrATVT i 

Fe(Ni) Si(Ti)Al 
^Co^Al Si(Al) 

Co(Ni)Ti 


Co(Ni,Al)Ti 

T?f 


__SKTj)_ 
■ Si(Aiyff 


Figure 13: Energy of formation (in eV/atom) for the Heusler alloys shown in Fig. 12 for Hf, Zr, 


For such a complex system, the energy spectrum introduced in previous sections becomes per- 
haps the only way to investigate the trends exhibited by these systems. The energy spectrum cor- 
responding to these numerical results is shown in Fig. 13. In spite of the similarities with the 
ternary case, it should be noted that the bonding scheme in the Heusler phase is quite different 
than that found in NiAl alloys. Fig. 13 shows that while Zr, Hf, Cu and Si show a distinct prefer- 
ence for the sublattice occupied by Ti and A1 atoms, there seems to be a low energy price to pay 
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for substitutions in A1 or Ti sites, even at the expense of an antistructure defect created when Ti 
goes to an A1 site in the X(Ti)Al cases. A clear distinction is seen for Co and Fe, where both show 
a marked preference for Ni sites. 

A complete presentation of the BFS results for all possible additions to all the Heusler alloys 
which can be obtained from the 16 elements considered in this work is, as mentioned above, 
beyond the scope of this paper, but not the computational method. We therefore limit our com- 
ments on such systems to the N^AlTi +X case (Table 7, Figs. 12-13) and to providing approxi- 
mate expressions using the BF rule (6) for the changes in lattice parameter for alloying additions 
in the dilute limit for general Heusler alloys (A 2 BC)!_ x X x , 

2^«1+|U: 07) 

a H 

with 

_ h(zh (18) 

Y H 

where y f = Bp ] , k, = , for i= A, B, C or X, and where y H and X H are given by 


1h = 


2Y a +Yb + Yc 
4 


(19) 


and 


x H .l bilhlh. ( 20 ) 

Several intermetallics of potential industrial interest form Heusler precipitates for small levels 
of additions. Some alloying additions, like Hf, Zr, Ta or Nb in NiAl, have a low solubility limit 
and result in the precipitation of Heusler particles. Such effect was studied in detail for the case of 
Ti in NiAl (10), where BFS was applied to the determination of the solubility limit of Ti in NiAl 
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and second phase formation. It is then important to examine the role of other additions and to 
determine, beyond their site preference in either phase (in the case of two phase alloys), their par- 
titioning behavior. 

The current formalism can be easily extended to provide similar information for the large fam- 
ily of alloys that exhibit the Heusler structure (N^AIX, with X= Nb, Ta, Hf, Zr, C 02 AIX, 
X=Nb,Ta,Hf,Zr and Ti; Cu 2 A1X, X= Hf, Zr; and Fe 2 AlX, X= V, Cr, Co, Ni). In addition, as 
noted for the case of binary systems, there are no restrictions in the formalism for the number of 
simultaneous alloying additions considered in the base alloy. Natural extensions of this work to 
high order additions in binary and ternary ordered intermetallics as well as other base systems 
(Ni 3 Al, Ni-base superalloys, etc. ) are also being performed. 

Conclusions 

The BFS method was used to determine the site preference scheme for a large number of 
alloying additions to NiAl, FeAl and CoAl. A detailed analysis of all the possible site substitu- 
tions was performed enhancing our understanding of the interrelationship between multiple alloy- 
ing additions. Because the agreement with experiment was excellent in all cases, the analysis was 
extended to ternary (Heusler) alloys where no data on the effect of alloying additions on structure 
exist. 
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